Fluctuation formula in the Nose-Hoover thermostated Lorentz gas 



M. DolowschialE and Z. KovacfQ 

Institute for Theoretical Physics, Eotvos University, Pf. 32, H-1518 Budapest, Hungary 

In this paper we examine numerically the Gallavotti-Cohen fluctuation formula for phase-space 
contraction rate and entropy production rate fluctuations in the Nose-Hoover thermostated periodic 
Lorentz gas. Our results indicate that while the phase-space contraction rate fluctuations violate 
the fluctuation formula near equilibrium states, the entropy production rate fluctuations obey this 
formula near and far from equilibrium states as well. 



PACS numbers; 05.45.-a, 05.70.Ln 



I. INTRODUCTION 

In recent years a large number of papers focused on 
various fluctuation formulas (FFs) with both theoreti- 
cal and numerical tools, and today these seem to be one 
of the most interesting results in the field of statistical 
physics of nonequilibrium systems P| . This behavior was 
observed numerically in a system of thermostated fluid 
particles undergoing shear flow (21 . It was an important 
property of the FF that it seemed to be valid for large 
external forcings as well 's^, therefore it was considered 
that it could shed some light on the thermodynamical 
behavior of systems far from equilibrium. Consequently, 
significant theoretical efforts have been made to find a 
common property behind the observed FF and these ef- 
forts resulted in various fluctuation theorems (FTs): the 
Gallavotti-Cohen approach built on the chaotic hypoth- 
esis 0,|3jthe Evans-Searles theorems @, deterministic 
local FT (3, the FT for stochastic systems 1l9l[T?]|. the 
theorem of Maes established on the Gibbs Property , 
and the FT for open systems The Gallavotti-Cohen 
FT/FF serves as the basis of the numerical investigations 
presented in this paper. 

All of the applied theoretical methods shared the prop- 
erty of putting extra presumptions on the physical sys- 
tems (i.e., chaoticity, stochasticity) that could not be 
proved a priori. This situation naturally raised the need 
to study the FF numerically and compare the numerical 
results to the theoretical predictions. Up to now several 
physical models have been investigated numerically, such 
as the two-dimensional (2D) reversibly damped fluids 7], 
the chains of weakly interacting cat maps 14j , the Fermi- 
Ulam-Pasta chain 13], and the periodic Lorentz Gas 
(PLC) thermostated by the Gaussian isokinetic (GIK) 
thermostat HI HI. 



The FF is a symmetry property of the probability den- 
sity function (PDF) of a dynamically measured quan- 
tity TT connecting the probabilities of measuring tt values 
with equal magnitudes but opposite signs. More pre- 
cisely, let T^rit) denote the quantity p averaged over a 



time interval of length r centered around time t: TTrit) — 

r/2 

^ / p{t + t')dt'. Considering it as a stochastic vari- 

-r/2 

able X, its statistical properties in a steady state can be 
characterized by the PDF nT-(a;). The FF states that the 
PDF IIt{x) has the following property: 



hmimil^ 



(1) 



In other words, for large enough r values the probabil- 
ity of observing — tt,- is exponentially smaller than the 
probability of observing tTt- . 

In p^ . Il6| | the examined physical quantity tt^. (t) was 
the phase-space contraction rate (PSCR), which due to 
this special property of the applied GIK thermostat was 
equal to the thermodynamical entropy production rate 
(EPR) at any given time. However, this identity does 
not hold in general 0|, and PSCR and EPR fluctua- 
tions can have different PDFs as, e.g., in the case of 
Nose-Hoover thermostated systems (see Sec. EJ. The 
theoretical methods using the Sinai-Ruelle-Bowen (SRB) 
measure to calculate the probabilities of trajectory 
segments predict that the quantity obeying the FF is 
the PSCR and further nontrivial theoretical efforts are 
needed to establish similar results for the EPR fluctua- 
tions p7j . 

The main purpose of examining the Nose-Hoover (NH) 
thermostated PLC is to investigate which of the above- 
mentioned physical quantities obey the FF in a system, 
where the corresponding PDFs are not identical. In addi- 
tion to this, checking the NH thermostated PLC numer- 
ically against the FF is in itself an important task, given 
that the NH thermostat is one of the two generally used 
dynamical thermostats; the transport properties of this 
model have been recently investigated in 18]. In Sec. lU 
we describe the examined model; in Sec. IIIII we present 
our numerical results; and in Sec. II VI we summarize our 
conclusions. 



II. THE SYSTEM 
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One of the most investigated models suitable for study- 
ing transport phenomena is the field-driven thermostated 
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periodic Lorentz Gas. This model consists of a charged 
particle subjected to external electric field and moving 
in the lattice of elastic scatterers. Due to the applied 
electric field, one must use a thermostating mechanism 
to achieve a steady state in the system; such a tool is 
a dynamical thermostat. Two types of dynamical ther- 
mostats have been applied to the PLG up to now: the 
GIK thermostat producing microcanonical distribution 
[Toj and the Nose-Hoover thermostat producing canoni- 
cal distribution in equilibrium 0, |23| • 

We present the equations of motion of the two- 
dimensional NH thermostated PLG in dimensionless 
form: mass and electric charge are measured in units of 
the particle's mass m and charge q and the unit length 
is chosen to be the radius of the scatterers {R = 1). Let 
q = {qi, (72) denote the position and p = (pi,P2) the mo- 
mentum of the particle, then the phase-space vector of 
the system is F = (q, p, (), where ( is the state variable 
of the thermal reservoir. Between two subsequent colli- 
sions the state of the system is evolved smoothly by the 
differential equation 



'^r{t) - i/i C{t)dt' and = i Sil^dt 
and check the validity of the FF for them. In order to 
perform this task we should evolve the state of the system 
r along a long trajectory, which requires the algorithm to 
be very efficient. This need motivated us to implement an 
event driven algorithm that generates and handles events, 
such as the collision of the particle with a scatterer and 
the replacement of the particle from one simulation cell 
into the other. The most sensitive issue when applying 
such an algorithm is to determine the point of time when 
a specific event occurs; in computer science this problem 
is known as collision detection. Since in the case of the 
Nose- Hoover thermostat the velocity of the particle is not 
upper bounded, we could not have chosen the simplest 
such method, the so-called naive algorithm, which could 
have been applied in the case of the GIK thermostat. 
Instead of this we have applied the method of building the 
time estimation into the Runge-Kutta integrator, which 
is used in simulating particle laden flows and coupled 
particle-field systems (see Ref. [2l|'). 
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and is transformed abruptly at every elastic collision. In 
this equation E is the the external electric field, Tj-esp is 
the response time of the reservoir, and T is the temper- 
ature satisfying (p^) = 2T. 

In the simulations presented in this paper we have used 
a square lattice of circular scatterers, however, we have 
investigated numerically other lattices as well (e.g, trian- 
gular), but have not found any relevant differences con- 
cerning the results presented in Sec. IIIII 

Energy dissipation can be measured by the phase-space 
contraction rate cr and can be computed by taking the 
divergence of the right-hand side of Eq. (0) as: 



a{t) = -divf (t) = C{t) 



(3) 



The entropy production rate ^ can be formally defined 
by the expression of irreversible thermodynamics 
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J(t)E _ p(t)E 
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(4) 



where T is the kinetic temperature. We note that this 
quantity is identical to the dissipation function of the 
Evans-Searles theorem It can be shown that in this 
model (cr) = (^), however, the identity cr(t) = ^(i) does 
not hold at all times, as opposed to the case of the GIK 
thermostated PLG. 



III. NUMERICAL RESULTS 

The objective of the numerical simulation is to measure 
the PDFs Tit{x) and ^^(a;) of the averaged quantities 
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FIG. 1: The probability density function T,-r(x) of the aver- 
aged phase space contraction rate fluctuations o"t- in a con- 
figuration close to equilibrium (E = (0.1,0.2)). 
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FIG. 2: The probability density function "Eri^x) of the aver- 
aged entropy production rate fiuctuations in a configura- 
tion close to equilibrium (E = (0.1,0.2)). 
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The PDFs 'Sr{x) and S^(a;) were constructed by pe- 
riodically computing the quantity a^-it) and £,rit) along 
a long particle trajectory and making a histogram of the 
computed data. In building the histogram we have used 
overlapping and nonoverlapping windows techniques as 
well, but have not found any relevant differences between 
them concerning the results of this paper. Throughout 
the presented numerical experiments we have used the 
T = 1.0 and Tresp = 1.0 values, however we have tested 
several other configurations as well. We have simulated 
t — 10^ long particle trajectories resulting in approxi- 
mately 2 X 10^ collisions with the scatterers. 

Figures El El a-nd 01 show the functional form of the 
PDFs St (2;) and 'E.r{x) for small and large external fields. 
Examining the figures one can make the following inter- 
esting observations: 

1. For low |E| values (close to equilibrium) in Figs lll2l 

the PDF T,r{x) seems to be more symmetric than 
the PDF Er{x). 

2. For low |E| values in Fig. [H and Fig. El the PDF 

'E.r{x) seems strictly narrower than ST-(a;); indeed 
in the |E| ^0 limit Er{x) should converge to a 
Dirac 6 function . [^t(0 = 0], opposed to dr that 
can fluctuate even in equilibrium. 

3. For high r values in Fig. El El and El the curves 
seem to be indistinguishable from a Gaussian; in- 
deed fitting a Gaussian onto the measured values 
yields an excellent visual agreement. 




FIG. 4: The probability density function S,t{x) of the aver- 
aged entropy production rate fluctuations in a configura- 
tion far from equilibrium (E = (0.8, 1.6)). 
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FIG. 5: The quantity D%(x) for the phase space contraction 
rate in a configuration close to equilibrium (E = (0.1, 0.2)). 



FIG. 3: The probability density function T,t{x) of the aver- 
aged phase space contraction rate fiuctuations CTt in a con- 
figuration far from equilibrium (E — (0.8, 1.6))]. 

With the measured values of [x) and (x) one can 
check the FF in this model. In order to visualize the FF 
we may introduce the quantities 



T Y.r{-x)' 

1 , ^r{x) 
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FIG. 6: The quantity D\{x) for the entropy production rate 
in a configuration close to equilibrium (E = (0.1,0.2)). 



which are shown on Fig. El El 13 and El With these 
quantities Eq. Q reads as \ivar-*oo D'^{x) = x and 
limi-^oo D\{x) = x. Examining Fig. El El El and El one 

can conclude that 

1. For small external forcings and for numerically 
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2. The fluctuations seem to obey the FF both for 
small and large external forcings. 

3. As the external forcing grows, the a-r fluctuations 
seem to obey the FF for large t values (Fig. O 
similarly to fluctuations. 

We note that we have examined several other config- 
urations and have found no signiflcant qualitative differ- 
ences in the observed behavior of the PDFs. 



IV. CONCLUSION 



FIG. 7: The quantity D%(x') for the phase space contraction 
rate in a configuration far from equilibrium (E = (0.8, 1.6)). 
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FIG. 8: The quantity D^{x) for the entropy production rate 
in a configuration far from equilibrium (E = (0.8, 1.6)). 



In this paper we have presented numerical evidence 
showing that the phase-space contraction rate fluctua- 
tions violate the fluctuation formula [Eq. |^] in or close 
to equilibrium in the Nose-Hoover thermostated periodic 
Lorentz gas. This observation is completely in line with 
the theoretical predictions of Evans et al in Ref. . 

On the other hand, we also demonstrated that the en- 
tropy production rate fluctuations satisfy the fluctuation 
formula, which is by no means trivial. It should also be 
noted that from the physical point of view the entropy 
production rate is the more relevant quantity due to its 
relation to thermodynamics and the availability for mea- 
surement in physical experiments. A direct consequence 
of our results is that the phase-space contraction rate 
and entropy production rate cannot be treated as inter- 
changeable qunatities in fluctuation formulas, therefore 
statements regarding the applicability of the FF for vari- 
ous models should always clarify which fluctuations they 
refer to. 



available r values the fluctuations seem to vi- 
olate the FF (Fig. El; this observation is supported 
by the fact that any symmetric function substituted 
in Eq. ||SJ) yields zero, and due to the time-reversing 
symmetry of Eq. ^ the PDF E^- {x) is expected to 
be close to a symmetric function for small external 
forcings (see Ref. 22]). 
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